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Abstract 

We show that the fohowing three systems related to various hydrodynamical 
approximations: the Korteweg-de Vries equation, the Camassa-Holm equation, 
and the Hunter-Saxton equation, have the same symmetry group and similar bi- 
hamiltonian structures. It turns out that their configuration space is the Virasoro 
group and all three dynamical systems can be regarded as equations of the geodesic 
flow associated to different right-invariant metrics on this group or on appropri- 
ate homogeneous spaces. In particular, we describe how Arnold's approach to the 
Euler equations as geodesic flows of one-sided invariant metrics extends from Lie 
groups to homogeneous spaces. 

We also show that the above three cases describe all generic bihamiltonian 
systems which are related to the Virasoro group and can be integrated by the 
translation argument principle: they correspond precisely to the three different 
types of generic Virasoro orbits. Finally, we discuss interrelation between the 
above metrics and Kahler structures on Virasoro orbits as well as open questions 
regarding integrable systems corresponding to a flner classiflcation of the orbits. 



1 Introduction 



One of the main mechanisms of integrability of evolution equations is the presence of two 
compatible Hamiltonian structures. In this paper we compare Hamiltonian properties of 
three extensively studied nonlinear equations of mathematical physics, related to various 
hydrodynamical approximations: the Korteweg-de Vries equation 

Ut = -3uUx + CUxxx, (1-1) 
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the Camassa-Holm equation 



^txx 3uUx ~\~ '^UxUxx ~l~ ^^xxx ~l~ ^^xxx (-^■^) 



derived as a shallow water equation in |PH|| (see also the paper |[FF|| ), and the Hunter- 
Saxon equation | |H|j| 

'^txx '^^x^xx ^'^xxxy (-^■'^) 

describing weakly nonlinear unidirectional waves. All three equations are known to be 
bihamiltonian and to possess infinitely many conserved quantities, as well as remarkable 
soliton or soliton-like solutions. A motivation for our paper was the paper [ BSS ], which 
described scattering theory for all three equations in a unified way. 

As we show in this paper, the main reason why such a common treatment is possible 
is that all these equations have the same symmetry group. It turns out to be the 
Virasoro group, a one-dimensional extension of the group of smooth transformations of 
the circle. More precisely, the Virasoro group serves as the configuration space, and all 
three equations can be regarded as equations of the geodesic fiow related to different 
right-invariant metrics on this group (in the case of KdV and CH) or on an associated 
homogeneous space (in the case of HS). (Here we will be mostly concerned with the 
periodic case, though many statements can be extended to the case of rapidly decaying 
potentials on the real line.) 

The main goal of this paper is to give a description of the three equations as bi- 
hamiltonian systems on the dual to the Virasoro algebra, and to relate them to the 
geometry of the Virasoro coadjoint orbits. One of the corresponding Hamiltonian (or 
Poisson) structures is provided by the linear Lie-Poisson bracket, and it is the same for 
all three equations. The other Poisson structure is constant and can be viewed as a 
"linear structure frozen at a point." The corresponding "freezing points" are different 
for each equation. We will also see that, in a sense. Equations (|1.1|) , (|1.2| ), and ( |T73|) 
exhaust all generic possibilities, and among them the Camassa-Holm equation ( |1.2| ) is 
the "most general" equation that can be obtained by the "freezing argument" method 
on the dual Virasoro space. 

We tried to make the paper self-contained, including in it necessary background on 
the Euler equations and the classification of Virasoro orbits. For additional information 
we refer the interested reader to the expositions in |^eg|| or ||AK|| , as well as to the original 
papers listed in the bibliography. 

In ||Arn|| V. Arnold suggested a general framework for the Euler equations on an 
arbitrary (possibly infinite-dimensional) group, which we recall below. In this framework 
the Euler equation describes a geodesic fiow with respect to a suitable one-sided invariant 
Riemannian metric on the given group. 

In Section 2 we show how Arnold's approach to the Euler equation works for the 
Virasoro group and provides a natural geometric setting for the Korteweg-de Vries and 
Camassa-Holm equations. In Section 3 we give a Hamiltonian reformulation of the Euler 
equation. In Section 4 we extend this approach to include geodesic fiows on homogeneous 
spaces and then use it to describe the Hunter-Saxton equation and its relatives.]^ This 



^In particular, the Harry Dym equation [Kru|, |HZb] can be found as one of the equations in the 
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extension might be thought of as a version of the Hamiltonian formahsm for homogeneous 
spaces developed in |[rh|| , |pS|| , which is apphed to the case of a degenerate metric and 



in our infinite-dimensional situation. 

In Sections 5-6 we develop the bihamiltonian formalism for the Euler systems. We 
show why the three equations above represent three main classes (and exhaust the 
rotation-invariant) bihamiltonian systems on the Virasoro algebra that can be integrated 
by means of the "freezing" (called also, the translation of) argument method. It turns out 
that the above three equations correspond to three different types of Virasoro coadjoint 
orbits of low codimensions. 

In our classification of Poisson pairs and the corresponding equations we rely heavily 
on the classical classification of the Virasoro orbits, and we recall it in the Appendix. 
An interesting open question is to extend the classification of the equations to orbits 
of higher codimension, as well as to show how the discrete invariant of Virasoro orbits 
manifests itself in the related bihamiltonian systems. 



2 The Euler properties of the KdV, CH, and HS 
equations 

The main objects in our consideration will be the group Diff(S'^) of all diffeomorphisms 
of a circle, its Lie algebra vect{S^) of vector fields, and their central extensions. The 
following nontrivial one-dimensional extension of the algebra of vector fields has a special 
name. 

Definition 2.1 The Virasoro algebra is an extension of the Lie algebra vect(S'"'^) of 
vector fields on the circle: 

vir = vect{S^) © M 

with the commutator between pairs (consisting of a vector field and a real number) given 
by 

[{v{x)d^,b),{w{x)dx,c)] = ( {-vwcc + Vxw){x)d^, C{vdx,wdx) ) 



where 



C{vda:,wdx) = / VWxxxdx 



is the Gelfand-Fuchs cocycle. (Here x is a coordinate on the circle, the subscript x stands 
for the derivative in x, and denotes the vector field on S^.) 

Define the following two-parameter family of quadratic forms, "if^ ^-energies," on 
the Lie algebra vir: 

{{v{x)dx,b),{w{x)dx,c))H^^ = / {avw + l3vxWx)dx + bc. (2.4) 



bihamiltonian hierarchy associated with the HS system. This equation was also considered in | BSS| . In 
this way, the Harry Dym equation also becomes associated to the geodesic interpretation. 
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The case a = 1, /? = corresponds to the inner product, while a = (3 = 1 corresponds 
to H^. Given a 7^ and any /5, extend the if^ ^-energy to a right-invariant metric on 
the Virasoro group Vir. This group corresponds to the Virasoro algebra and is defined 
as follows. 

Definition 2.2 The Virasoro group Vir is the product 

Vir = Diff(5^) x M, 
where the group multiplication between the pairs is given by 

{^(x), a) o (0(x), 6) = ( o 0)(x), a + b + B{i;, 0) ) 

and 

B{tlj,(f)) := / log((^o0ydlog(0j 

is the Bott cocycle. (Here the diffeomorphisms of are described by functions, e.g., 
X (— > (p{x).) 

Having equipped the Virasoro group with those right-invariant metrics, one can con- 
sider the geodesic fiows they generate. 



Theorem 2.3 1) |PK|| The KdV equation is the Euler equation, describing the geodesic 
flow on the Virasoro group with respect to the right-invariant Li^-metric. 

2) ||M| The CH equation is the Euler equation for the geodesic flow on the same group 



with respect to the right-invariant Sobolev -metric. 

It turns out that one can give a similar description of the Hunter-Saxton equation as 
a geodesic flow on a homogeneous space related to the Virasoro algebra. Consider the 
-quadratic form (which is the if^ ^ -form with a = 0, (3 = 1) on the Virasoro algebra: 



{{v{x)d^,b),{w{x)d:r,c))Hi = / v^w^dx + bc. (2.5) 

Although this form is degenerate, as is the corresponding right-invariant metric on the 
Virasoro group, one can define a nondegenerate metric by descending on an appropriate 
quotient space. 

Theorem 2.4 The HS equation is the equation describing the geodesic flow on the ho- 
mogeneous space Vir /'Roi{S^) of the Virasoro group modulo rotations with respect to the 
right-invariant homogeneous metric. 

Note that one can also obtain the HS equation by considering the smaller homoge- 
neous space Diff(S'^)/Rot(S'^) of all diffeomorphisms of the circle modulo rotations, as 
we explain below. 

These three equations essentially exhaust the list of integrable systems associated 
with the Virasoro algebra and integrated by the freezing argument method, as we discuss 
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below. Note that their degenerations include, e.g., the inviscid Burgers equation (it 

corresponds to the L^-metric on the "centerless Virasoro" group, Diff(5'^)). 

2.5 Remark. Before proving the theorems, we recall the general set-up for the Euler 



equation on an arbitrary Lie group, suggested by V. Arnold in [|Arn|| . Consider a (possibly 
infinite-dimensional) Lie group G, which can be thought of as the configuration space 
of some physical system. (Examples from ||Arn|] : 5*0(3) for a rigid body or the group 



SDiff(M) of volume-preserving diffeomorphisms for an ideal fluid filling a domain M.) 
The tangent space at the identity of the Lie group G is the corresponding Lie algebra 
g. Fix some (positive definite) quadratic form, the "energy," on q. We consider right 
translations of this quadratic form to the tangent space at any point of the group (the 
"translational symmetry" of the energy). This way the energy defines a right-invariant 
Riemannian metric on the group G. The geodesic flow on G with respect to this energy 
metric represents the extremals of the least action principle, i.e. the actual motions of 
our physical system.^] 

To describe a geodesic on the Lie group with an initial velocity v{0), we transport 
its velocity vector at any moment t to the identity of the group (by using the right 
translation). This way we obtain the evolution law for v{t), given by a (non-linear) 
dynamical system dv/dt = F{v) on the Lie algebra q. 

Definition 2.6 The system on the Lie algebra g, describing the evolution of the velocity 
vector along a geodesic in a right-invariant metric on the Lie group G, is called the Euler 
equation corresponding to this metric on G. 




Figure 1: The vector v in the Lie algebra g traces the evolution of the velocity vector of a geodesic 7 
on the group. The inertia operator A sends w to a vector m in the dual space g*. 

In particular, the above scheme works for the Virasoro group (see Theorem p^.3| ) and 
allows one to describe the Korteweg-de Vries and Camassa-Holm equations as geodesic 
equations on that group. It also can be extended to include geodesic flows on homoge- 
neous spaces and to describe the Hunter-Saxton equation, as we discuss below. 



^For a rigid body one has to consider left translations, but in our exposition we stick to the right- 
invariant case in view of its applications to the groups of diffeomorphisms. 
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3 Hamiltonian framework for the Euler equations 



We start with preliminaries on Lie algebras and Poisson structures. 

Definition 3.1 The dual space g* to any Lie algebra g carries a natural Lie-Poisson 
structure; 

{f,9}Lpim) := {[df,dg],m) 

for any m G g* and any two smooth functions f,g on g*. (Here the differentials are 
taken at the point m, and (-, ■) is a natural pairing between g and q* .) 

In other words, the Lie-Poisson bracket of two linear functions on g* is equal to their 
commutator as elements of the Lie algebra g itself. 

Proposition 3.2 The Hamiltonian vector field on g* corresponding to a Hamiltonian 
function f and computed with respect to the Lie-Poisson structure has the following form: 

dm 

— = ad^fm. (3.6) 

Proof. Let dm/dt = Xf be the corresponding Hamiltonian field. Then for any 
function g G C°°(g*) one has the identities 

ixfdg\^ = Lxj,g\m = {f,9}Lpim) = {[df,dg],m) = {dg,&d*ajm). 

This implies that Xf = ad^^m. □ 

3.3 Remark. The differential-geometric description of the Euler equation as a geodesic 
flow on a Lie group has a Hamiltonian reformulation. 

Fix the notation E{v) = ^{v, Av) for the energy quadratic form on g which we used 
to define the Riemannian metric. Identify the Lie algebra and its dual with the help of 
this quadratic form. This identification A : q ^ g* (called the inertia operator) allows 
one to rewrite the Euler equation on the dual space g*, see Fig.l. 

It turns out that the Euler equation on g* is Hamiltonian with respect to the Lie- 
Poisson structure [ |Arn|] . Moreover, the corresponding Hamiltonian function is minus 
the energy quadratic form lifted from the Lie algebra to its dual space by the same 
identification: —H{m) = —^{A~^m,m), where m = Av. Here we are going to take it as 
the definition of the Euler equation (we use the Proposition above and the observation 
dH{m) = A'^m). 

Definition 3.4 The Euler equation on g* corresponding to the Hamiltonian —H(m) = 
— ^{A^^m,m) is given by the following explicit formula: 

dm _ 

as an evolution of a point m G g* . 
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3.5 Remark. The underlying reason for the Hamiltonian reformulation is the fact 
that any geodesic problem in Riemannian geometry can be described in terms of sym- 
plectic geometry. Geodesies on M are extremals of a quadratic Lagrangian (metric) on 
TM. They can also be described by the Hamiltonian flow on T*M for the quadratic 
Hamiltonian function obtained from the Lagrangian via the Legendre transform. 

If the manifold is a group G with a right-invariant metric then there exists the 
group action on the tangent bundle TG, as well as on the cotangent bundle T*G. By 
taking the quotient with respect to the group action, we obtain from the (symplectic) 
cotangent bundle T*G the Lie-Poisson structure on the cotangent space T*G\e = 0*, 
i.e., on the dual to the Lie algebra. The Hamiltonian function on T*G is dual to the 
Riemannian metric (viewed as a form on TG), and its restriction to g* is the quadratic 
form H{m) = ^{A~^m,m) , m G 0*. 

The geodesies of a left-invariant metric on G correspond to the Hamiltonian function 
H{m), while those of a right-invariant metric correspond to —H{m). 



Now we are ready to prove Theorem 2.3 on the Eulerian nature of the KdV and CH 



equations in the following slightly more general setting. 



Theorem 3.6 (= |2.3f ) The Euler equation describing the geodesic flow on the Virasoro 
group with respect to the right-invariant H\ ^ -metric with a ^ has the form: 

a{vt + 3vv^) - I3{v^^t + '^v^v^^ + w^xx) - bv^^^ = 0, (3.7) 

bt = 0. 

3.7 Remark. By choosing a = 1, (3 = one obtains the KdV equation, related to the 
L^-metric on the Virasoro algebra [PK|| . Similarly, for a = (3 = 1 one recovers a general 
form of the CH equation ||Mi|| . Note that by shifting v ^ f+const we get another form 
of the CH equation, in which the term Vxxx is replaced by Vx- Finally, if a = 0, /? = 1 
then Equation ( p. 71 ) becomes the HS (Hunter-Saxton) equation, which we discuss in the 
next section. 

The case 6 = corresponds to considering the non-extended Lie algebra vect{S^) of 
vector fields rather than the Virasoro algebra vir. Depending on the values of a and /3 
one obtains either the inviscid Burgers (also called, Hopf ) equation Vt + 'ivVx = or the 
non-extended CH equation. 

Proof of Theorem |3.6| . Recall that the Virasoro coadjoint action can be computed 
as follows. Let {{u{dx)'^, a)| m G C°°(S'), a G M} be the dual space to the Virasoro 
algebra with the natural pairing given by 

{{u{dxY, a), {wdxi c)) = I uw dx + ac. 



(Here we denote by u{x){dxY (or, shorter, by u{dxY) a quadratic differential on the 
circle.) The coadjoint operator is defined by the identity 

(ad(^9.,fe)(«(c^a;)^ a), {wdx,c)) = {{u{dxf , a), [{vdx,b), {wdx,c)]). 
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Using the definition of tlie Virasoro commutator and integrating by parts we obtain that 
the right- hand- side is equal to 



w{2uvx + UxV — aVxxx) dx. 

1 

Thus the coadjoint operator is 

ad*^Q^ ,,){u{dxf, a) = i{2uVx + UxV - aVxxx){dxf, 0). (3. 

The Hl^ p -energy ( p.4|) on the Virasoro algebra 



{{vdx, b), {wdx, c)) = / {avw + (3 VxWx) dx + be = v Aw dx + be. 

corresponds to the general inertia operator A : vir — >■ vir*, given by 

{vdx,b)^{{Av){dxf, b), 

where A := a — is a second order differential operator. This operator is non- 
degenerate on vir for a ^ 0, while for a = it has a non-trivial kernel consisting of 
constant vector fields on S^. 
Now the Euler equation 

^{u{dxY,a) = -ad*^-i(u(dxr,a)iu{dxy,a) 



on vir* (see Definition ( |3.4| )) assumes the form 



— {u{dxY^a) = — ((2'uA ^Ux + UxA — aA ^ Uxxx){dx)'^ ., O) . 

(Here we substituted {vdx,b) = A^^{u{dxY,a) = {{A~^u){dxY,a) into the expression 
for ad*.) 

In terms of f = A^^m (the first component of) this equation becomes 

A(At>) = -2{Av)vx - iAvx)v + bvxxx, 
at 

which is equivalent to the equation (p.7|), since A = a — pd^. For the second component 
we find that b does not change in time: bt = 0. □ 

3.8 Remark. In the proof we assumed that the inertia operator is invertible. In the next 
section we discuss the precise relation of the geodesic and Hamiltonian approaches in the 
case of a degenerate metric and show what reductions are necessary for the corresponding 
Euler equation to make sense. 
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4 The Euler equations on homogeneous spaces 



Let G be a Lie group and K its subgroup. Consider the space GjK of right cosets 
{Kg I g G G}. Then the group G acts on them on the right. Here we are going to 
develop the formahsm for the Euler equation, describing the geodesic flow on G/K with 
respect to a right-invariant metric. 

One immediately encounters the following difficulty: not every right-invariant metric 
on the group G, degenerate along K at the identity, descends to a metric on the space 
of right cosets G/K (see Example ([4.3| a) below). To formulate the condition, which the 
degenerate metric should satisfy, let us consider the corresponding problem at the level 
of Lie algebras. 

Let g be a Lie algebra, and A : g q* a degenerate inertia operator. Suppose that 
the kernel of A is a Lie subalgebra 6. (In other words, the corresponding energy form 
E{v) = ^{v,Av) vanishes for all v G t.) Consider the right-invariant degenerate metric 
Eq on the group G obtained by translating the quadratic form E from identity to any 
point of the group. 

Theorem 4.1 The right-invariant form Eg on a group G descends to a form on the 
space G/K of right cosets if and only if the quadratic form E on the Lie algebra g 
vanishes on the Lie subalgebra 6 and is Ad-invariant with respect to the action of this 
subalgebra. 

4.2 Remark. The condition of Ad-invariance for E reads as follows: 

{adwV,Au) = — {v , A{a.dwu)) 

for all w G t and any u,v G Q. The above is an infinitesimal version of the invar iance of 
E with respect to the subgroup action: 

{{Adkv),A{Adku)) = {v,Au) 

for all k e K. 

4.3 Example, a) (Rotations of a rod.) The configuration space of a rod in fixed at 
its center of mass is S*^. It can be obtained from the configuration space of a rigid body 
by moding out rotations about one of its axes: S'^ = \ 5*0(3). 

Suppose that A : M'^ is a degenerate inertia operator with one vanishing eigen- 

value. The corresponding eigenvector generates the 1-dimensional rotation subgroup 
in S0{3). It is not difficult to see that the bi-invariance condition imposes the following 
restriction on A: its non-vanishing eigenvalues must be equal. (Indeed, the S'^-action 
sends one of these two eigenvectors to the other.) Then the corresponding degenerate 
metric on 5*0(3) descends to (a multiple of) the standard metric on the sphere, which 
is the space of cosets: 5^ = 5^ \ 50(3). Geodesies with respect to the standard metric 
on 5*^ are the great circles. These geodesies describe all free motions of the rod.0 (We 

^ For the rod, as well as for a rigid body mentioned above, we consider the left-invariant metrics, 
and hence, left cosets \ 50(3). For the group of diffeomorphisms we study the right-invariant metrics 
and the space of right cosets Diff (5^)/Rot(5'^). 
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can see that the only parameter of the rod is its length. In terms of the inertia operator, 
this corresponds to the choice of the nonzero eigenvalue.) 

Note that the inertia operator A = diag(A,/i, 0) with A 7^ /i does not correspond to 
any physical object. The corresponding degenerate metric on S0{3) is not S'^-invariant, 
and hence it does not descend to S'^. 

b) (The Hunter-Saxton equation.) Consider the group of diffeomorphisms Diff(5'^) 
and its quotient Diff (S'^)/Rot(5'^) by the subgroup of rotations Rot(5'^). Consider the 
degenerate quadratic form on the corresponding Lie algebra vect{S^) given by the ho- 
mogeneous if^-energy: 

Eo{vdx) = 7: v^dx. 



51 



(Similarly, for the Virasoro algebra, consider the energy Eq^vOx, h) = ^(/^i dx + IP').) 

This energy vanishes on constant vector fields. Those fields generate the subgroup 
Rot(S'^) of rotations of the circle S^. One can see that the form EQ{vdx) is bi-invariant 
with respect to the circle action, since the energy is invariant with respect to translations 
X I— > x+ const. (The same holds for the Gelfand-Fuchs cocycle, and the energy Eq^vOx, h) 
on the extended algebra.) Hence the corresponding right-invariant metric on Diff(S'^) 
descends to the quotient Diff(5'^)/Rot(5'^). We will see that the geodesies with respect 
to this metric are described by the HS equation. 

4.4 Remark. One hopes that a certain modification of this approach can be applied 
to generalized flows in [[Bre|| , where the fluid particles in 3D can move freely and inde- 



pendently along one coordinate. The corresponding subgroup K here might be that of 
fiberwise diffeomorphisms along a coordinate. 



Proof of Theorem [4.1| . First of all we note that the quadratic form E owq induced 
from a nondegenerate form on the quotient g/t is degenerate exactly along 6 (i.e., 6 is 
its null subspace). 

Let 6 be a Lie subalgebra and consider the quotient of the corresponding groups 
G/K. Note that the restriction of the energy form to the subgroup K is zero. Indeed, 
the latter is nothing but the right translation from the identity of the energy form on 
the subalgebra 

We would like to compare the right-invariant metric Eq on G at two different points 
kiQ and k^g of the same coset Kg. Then the element k = k2k^^ G K sends ki to 
^2 by means of the left translation. This translation also identifies the tangent spaces 
to G along the same coset Kg, see Fig. 2. The energy form Eq is invariant under this 
identification, since it is bi-invariant with respect to the action of elements of K. Finally, 
note that the energy Eg is degenerate along cosets. (Indeed, it vanishes on K, the 
"identity coset," and it is invariant with respect to right translations, which shuffle the 
cosets.) 

Therefore, the corresponding energy form descends to the coset space. □ 

4.5 Remark. From the Hamiltonian point of view, the geodesic picture on a homo- 
geneous space corresponds to a Hamiltonian reduction of the non-degenerate case with 
respect to the subgroup K action. 
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Figure 2: Defining a right-invariant form on the space of right cosets G/K. 

More precisely, consider the cotangent bundle T*(G/K) of the metric space G/K. 
Similarly to the non-degenerate case, look at the fiber {Q/t)* over the "identity coset" 
K. This space (0/6)* can be naturally identified with the image L = imA C Q* of the 
degenerate inertia operator A : g — »• g*, for which 6 = keiA, (or, equivalently, with the 
subspace L = C g*, the annihilator of Hn g*). 

The subgroup K, being a stabilizer of the "identity coset" and hence of L = {g/t}*, 
acts on L respecting the Poisson structure. Therefore, there is a natural Poisson structure 
on the quotient L/K := imA/ad*^^. 

Furthermore, one can define the corresponding Hamiltonian function on L by the 
same formula as above: 

HL{m) = - {m,A-^m). (4.9) 

This function is ad*i^-invariant and hence is well-defined on the quotient L/K. 

Suppose A generates a right-invariant and ad /^-invariant metric on the group G (i.e., 
satisfies the conditions of Theorem [4.1| ), so that it makes sense to consider geodesies of 
the corresponding metric on G/K. 

Then the above consideration provides the following limiting (degenerate) case of 
Arnold's theorem (cf. Definition |3.4|) . 

Theorem 4.6 The Euler equation, which corresponds to the inertia operator A and 
describes the geodesic flow on the space G/K of right cosets, has the following Hamilto- 
nian form on L/K = imA/ad*K : it is the quotient with respect to the K-action of the 
restriction to L G Q* of the following Hamiltonian equation on g* 

— = -ad^-i^m 

for m E L = imA. 

Now we are ready to complete the argument showing the Eulerian nature of the 
Hunter-Saxton equation. 

Theorem 4.7 (=|2T^) The HS equation 
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is a well-defined equation on the equivalence classes of periodic functions 

{v{x) v{x + p) + q for any p, q}. 

It describes the geodesic flow on the homogeneous space Vir /KotlS^) of the Virasoro 
group modulo rotations with respect to the right-invariant homogeneous metric. 

Proof of Theorem |4.7| . The homogeneous if^-metric on the Virasoro algebra is 
related to the degenerate inertia operator A : vir —>■ vir* sending 

[vd.^b)^{-{dlv){dx)\b). 

Its image L C vir* consists of pairs {u{dxY, a), where functions u have zero mean: 

[u{dxY^a) I u{x) = 



The action of the subgroup K = Rot(S'^) identifies those functions u that differ by a 
rotation: u{x) ~ u{x + Thus we come to a Hamiltonian equation on L/K. 

For explicit calculations, recall that the present case corresponds to setting a = 0, 



/? = — 1 in the proof of Theorem |3.6| . Choosing these values in Equation ( p.7|) , one arrives 
at the HS equation 

on vir. In order to obtain equation describing the evolution of v rather than that of Vxx 
we observe that 

"^VccVxx + vv^xx + bv^^^ = (vv^x + ^{Vxf + bVxx 
and hence, integrating both sides in x, we obtain 

vtx = -vv:cx - - bv^^ + r, 

where r is an arbitrary constant. This constant is uniquely determined by the condition 
that the right-hand-side is a complete derivative, i.e., by ^^^{{vxY /2 + rQ)dx = 0. Then 

vt = -vv^ + d~^{{v^f/2 + ro) - bv^ + g, 

where q is an arbitrary constant. 

Whence Vt (and hence the evolution of v) is defined only up to addition of a pa- 
rameter q and a multiple of Vx. (The latter is the velocity of the rotation subgroup: 
Vx = ^|p=o'^(2; + p).) This manifests the fact that the evolution of v is defined on the 
equivalence classes {v{x-\-p) + q}. The inertia operator A sends this equation on classes 
to the Hamiltonian equation on the quotient L/K. 

Note that the equivalence classes above absorb the term bvx (respectively, bvxxx for 
Vtxx'- consider a shift v ^ v -\- const) and one obtains the HS equation in its standard 
form (O). □ 



One should notice that setting 6 = corresponds to the Euler equation on the quo- 
tient Diff (S'^)/Rot(S'^). Thus in the homogeneous case the consideration of the central 
extension does not give anything new, since the Euler equations for Diff (S'^)/Rot(S'^) 
and V"ir/Rot(5'^) are equivalent. 
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5 Bihamiltonian structures for the equations 



To formulate our next result we need to recall some generalities on bihamiltonian systems. 

Definition 5.1 Assume that a manifold M is equipped with two Poisson structures 
{., .}o and They are said to be compatible (or, form a Poisson pairj if all of 

their linear combinations {., .}o + A{., .}i are also Poisson structures. 

A dynamical system dm/dt = F{m) on M is called bi-Hamiltonian if the vector field 
F is Hamiltonian with respect to both structures {., .}o and {., .}i. 

Consider the dual space q* to a Lie algebra q. As we discussed above, it is equipped 
with the Lie-Poisson structure: 

{f,9}Lpim) := {[df,dg],m) 

where m G g* and f,g are two arbitrary functions on q*. 

Now fix a point mg in q*. One can associate to this point another Poisson bracket 
on 0* as follows, see ||Ma|| . 

Definition 5.2 The constant Poisson bracket associated to a point mo G g* is the 
bracket {., .}o on the dual space Q* defined by 

{/,^}oM := {W,d9],mo) 

for any two smooth functions f, g on the dual space, and any m G g* . The differentials 
df, dg of the functions /, g are taken, as above, at the point m and are regarded as 
elements of the Lie algebra itself. 

The constant bracket depends on the choice of the "freezing" point mo, while the 
Lie-Poisson bracket is defined by the Lie algebra structure only. Note that the brackets 
{..I-} LP and {.,.}o coincide at the point mo itself, and, moreover, the bivector defining 
the constant bracket {., .}o is the same at all points m. 

Proposition 5.3 The brackets {•,.}lp ond {•,.}o are compatible for every "freezing" 
point niQ. 

Proof. Indeed, any linear combination {., .}\ := {., .j^p + A{., .}o is again a Poisson 
bracket, since it is just the linear Lie-Poisson structure {., .}lp translated from the origin 
to the point — Amg. □ 

5.4 Remark. Explicitly, the Hamiltonian equation on g* with the Hamiltonian function 
/ and computed with respect to the constant Poisson structure frozen at a point mo G g* 
has the following form: 

dm 

— = ad^^mo, (5.10) 



as a modification of Proposition |3.2| shows. 
Now we can formulate another main result. 
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Theorem 5.5 The Euler equation ^3. Tj ) for the H^p -metric (with a 0) on the Vira- 
soro group is hihamiltonian on the dual vir* of the Virasoro algebra. The corresponding 
"freezing" point in vir* is {^{dx)'^, (3). 

In Appendix we show that the dual space vir* = {{u{x){dx)'^ , a)} can be thought of 
as the space of HiU's operators {—ad^ + u{x)}. In these terms the above theorem can 
be stated as follows: The ^ -metric on vir given by the inertia operator A = a — (3dl 
is hihamiltonian on vir* with "freezing" at the point a/2 — fSd^. 

5.6 Remark. The KdV and CH equations are hihamiltonian on the Virasoro dual. The 
corresponding "freezing" points (uo((ia;)^, Oq) in vir* are {uq = 1/2, = 1) for the CH 
equation and (mq = 1/2, = 0) for the KdV equation, see Fig.3, as they are related to 
the H^- and L^-energies, respectively. 

To describe hihamiltonian nature of the HS equation on the reduced space L/K, 



discussed in Theorems [4.6| - [4.7| , one should consider the following analog of the con- 
stant Poisson structure. Take the Lie-Poisson structure on vir* "frozen" at the point 
(mq = 0, oq = 1) and then push it forward to the corresponding quotient space for the HS 
equation. An alternative way to show integrability (rather than the hihamiltonian prop- 
erty) of this equation is to use integrability of CH and an infinite-dimensional version of 



the formalism developed in ||l'h| | , [ |Clj| ] 



5.7 Question. Which metrics on the Virasoro group (or which quadratic forms on 
the Virasoro algebra) correspond to the hihamiltonian system on vir* with "freezing" 
at a point {uo{x){dxy,ao) for non-constant mo(x)?0 For which uo{x) are these metrics 
positive definite (i.e., Riemannian rather than pseudo-Riemannian)? 




Figure 3: Locations of the "freezing" points for the KdV, CH, and HS equations in the Virasoro 
dual vir*. 

Proof of Theorem |5.5| . Let F{u, a) be a function on vir* and let {vdx, b) = 

{SF/Su,6F/6a) be a (variational) derivative of F at {u{dx)'^,a). Then the Hamiltonian 
equation with Hamiltonian function F computed with respect to the constant Poisson 
structure "frozen" at {uQ{dx)'^ , ao) has the form 

^^{u{dxf,a) = adl^Q^^^){uo{dx)'^,ao) = {i2uoVa: - aot;^^^)(c?x)^ O) . 



^After this paper was submitted, I. Zakharevich found a formula for the corresponding quadratic 
form on the Lie algebra in terms of Bloch solutions of the operator —d^ + u{x), [Za|. 
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(Here we used Remark ^.4| , the explicit form (|3.8|) of ad* for the Virasoro algebra, and 
the fact that Uq = const.) 

For the first component of this equation one has 

Setting Uq = a/2 and qq = P we obtain 2uo — aod^ = a — [3d1 = A, and this simplifies 
the equation to 

f-(^f))- 

To prove the Theorem one needs to show that for any a ^ and any f3 the Euler equation 
( ^.71 ) can also be expressed in the form ( p.ll| ) for an appropriate Hamiltonian function 
F. 

Next, consider the Hamiltonian function F of the form 



F{u, a) 



(The operator A := a — f3dl is invertible for a ^ 0.) By definition, the variational 
derivative {6F/6u,6F/6a) G vir of the functional F is determined by the following 
identity satisfied for any {C,{dx)'^,c) G vir*: 



f6F 




\ d 


V 5^' 


da J 


lie 



{i{dx)\ c), (—,— )) = — ^^F{u + ee, a + ec) 



Since we need only the partial variational derivative 5F /5u, we compute: 

F{u + ee, a) = 



d 
de 

d 



t=0 



de 



e=0 



J (^^{A~\u + eOr + \{A-\u + eOy{u + eO + ^{A-\u + e0.r'^ dx 
- j i - (^A-\A-\f + \{A-^uf + ]^k-\{k-^u)u) - aK-^u,}j dx. 



Thus, we have found that 



^ = - {aK-\k-\f + \{^-'uf + ]^k-\{k-^u)u) - aA-V. ) 



Now, we substitute the variational derivative 6F/6u into Equation (|5.11|) and then 
rewrite the obtained equation on the algebra vir, rather than on its dual vir*. The 
latter corresponds to rewriting the equation in terms of {vdx,b) = A~^{u{dx)'^,a), i.e., 
in terms of an unknown function v = A~^u and setting b = a. 
Finally, applying A to the equation we obtain 

A f ^ j = - (2avvx + ^A{vvx) + ^vAv^ + ^t;^At; - bv^:: 



15 



Recalling that A = a — [3d1 and collecting the terms we recover Equation (p.7|) . □ 

To explain in what sense the above equations are generic bihamiltonian systems on 
vir* (obtained by the freezing argument method), we need to consider symplectic leaves 
of the above Poisson structures. 



6 Hierarchies of Hamiltonians from compatible 
structures 

Recall that the symplectic leaves, i.e., maximal non-degenerate submanifolds, of the Lie- 
Poisson structure are the coadjoint orbits of the group action on q* (see, e.g., ||Kir|| or 
this also follows from Proposition |3.2| ). Therefore the functions constant on symplectic 
leaves (called Casimir functions) of the Lie-Poisson bracket are those functions on the 
dual space g* that are invariant under the coadjoint action. The tangent plane to the 
group coadjoint orbit at the point mo, as well as all the planes in q* parallel to this 
tangent plane are the symplectic leaves of the constant bracket frozen at the point mo . 

Definition 6.1 The codimension of the coadjoint orbit passing through rriQ will be called 
the codimension of the Poisson pair {., .}o and {., .}lp- 

It turns out that there are no Poisson pairs of codimension or 1 in the (smooth) 
Virasoro dual vir*, and that the Poisson pairs of codimension 2 can all be classified. 

Theorem 6.2 All Poisson pairs {., .}o (ind {., .}lp on vir* of codimension 2 belong to 
one of three classes according to the orbit type of the 'freezing" point {uo{dxy, ao). These 
classes can be represented by the points a) {{dx)^/2,l), b) (((ia;)^/2, 0), and c) (0,1). 

Proof. First we observe that for any Lie algebra g the list of coadjoint orbits in q* 
provides the list of normal forms for the constant and Lie-Poisson pairs as well. Indeed, 
let a point mo G g* be a normal form for all points mo belonging to the same coadjoint 
orbit as mo- This means that mo can be mapped by the group coadjoint action to its 
normal form: mo = Ad*mo. The group action Ad* : g* — g* is a linear operator on g*, 
which preserves the Lie-Poisson bracket on g*. It also maps the constant bracket frozen 
at mo to that frozen at mo. Thus the group action Ad* sends one Poisson pair to the 
other. 



The proof of Theorem |6.2| is based on the Virasoro orbit classification, which we recall 
in Appendix. Notice that a "cocentral value" ao is invariant on the orbits in vir* . This 
allows us to fix ao and consider the orbits in the hyperplane {{u{x){dxY , a)\ a = ao}. 

As shown in Corollary of Appendix, for ao 7^ there are exactly two types of 
Virasoro orbits of codimension 1 in this hyperplane that correspond to either a generic 
or a Jordan 2x2 block holonomy for Hill's operator — aoi9^ + u{x). If we discard the 
discrete invariant of these orbits (see Appendix), representatives for their normal forms 
can be chosen as stated in the a) and b) parts of the Theorem. 
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If Oq = there is just one orbit type of codimension 1 in the corresponding hyperplane, 
and represented by c), see Remark |7^. Note that in the whole dual space those orbits 
have codimension 2, taking extra dimension a into account. □ 

Given a Poisson pair one can generate a bihamiltonian dynamical system, by produc- 
ing a sequence of Hamiltonians in involution, according to the following Lenard scheme. 
Let {., .}a := {., ■}o + A{., .}i be the Poisson bracket on a manifold M for any A. Denote 
by hx its Casimir function on M parameterized by A. This means that {hx, f}x = for 
any function /. Expand hx in a power series: hx = h^ + hiX + where each coefficient 
hj is a function on M. The following theorem is well-known. 

Theorem 6.3 The functions hj, j = 1,2 . . . are Hamiltonians of a hierarchy of bihamil- 
tonian systems. In other words, each function hj generates the Hamiltonian field Xj with 
respect to the Poisson bracket {., .}i (i.e., Xj satisfies Lx^f = {hj, /}i for any f), which 
is also Hamiltonian for the other bracket {., .}o with Hamiltonian function — /ij+i (i.e., 
Lx f = —{hj+i,f}o for any f). Other functions hi, i ^ j are first integrals of the 
corresponding dynamical systems Xj. 

In other words, the functions hj, j = 0,1 . . . are in involution with respect to each 
of the two Poisson brackets {., .}o and {., 

Proof. Substituting the power series for hx into the Casimir condition we obtain: 

= {hx. f}x = {ho + h^X + /}o + X{ho + h^X + 

Collecting the terms at A°, A^, A^, . . . we obtain a sequence of identities: 

{ho,f}o = 0, {h^,f}o + {ho,f}i = 0, {h2,f}o + {hi,f}i = 0,... 

for any function /. The first identity expresses the fact that h^ is a Casimir function for 
the bracket {., .}o- The next one says that the Hamiltonian field for hi with respect to 
{., .}o coincides with the Hamiltonian field for — /iq and the bracket {•, - ji, and so on. 

To see that every function hi is a first integral for the equation generated by hj with 
respect to each bracket, we check that {hi, hj}k = 0, k = 0,1. Indeed, e.g., if i < j 

{hi, hj}i = — {/ij+i, hj}Q = {/ij+i, hj_i}i = . . . = 0, 

since we finally obtain the bracket (either {., .}o or {., of one of the functions hi with 
itself. □ 

Thus the choice of Casimir functions hx determines the corresponding (hierarchy of) 



dynamical systems. By combining Theorems |6.2| and |6.3| we get the following 



Corollary 6.4 The three types of Poisson pairs of codimension 2 on the Virasoro alge- 
bra correspond to the three integrable systems: CH, KdV, and HS. These three systems 
represent all generic Hamiltonian systems on vir* (modulo the ambiguity in the choice 
of Casimir), which can be integrated by the freezing argument method. 
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6.5 Remark. The corresponding "freezing" points in vir* represent all three types of 
the Virasoro coadjoint orbits of codimension 2. 

If the "freezing" point {uo^dx)"^ , ao) is generic, one obtains an equation "equivalent" 
to the CH equation. In this sense, the CH equation is the most general equation, which 
is encountered by applying the freezing argument method of integration; this is the case 
a) in Theorem |6.2i 

Two other equations can be recovered by confining the "freezing" point to special 
hypersurfaces in vir*. (In turn, these hypersurfaces are foliated into coadjoint orbits. 
Those orbits are of codimension 1 in the hypersurfaces, and hence of total codimension 
2 in vir*. The classification of the orbits will be discussed in detail in Appendix.) A 
generic point on the hyperplane qq = produces the KdV equation; see the case b) 
in Theorem |6.2| . The case c) in the same theorem corresponds to the HS equation if 
we consider the cone-like Virasoro orbits in the oq 7^ 0-case (see Appendix and Fig. 4). 
The latter (e.g., "freezing" at the point {uQ{dxY,ao) = (0, 1)) corresponds to the Euler 
equation with a degenerate metric on the group, which we discussed in Section |^. 

One could also consider a more subtle Virasoro orbit classification, where one dis- 
tinguishes between the two types of generic orbits in vir*: hyperbolic and elliptic ones, 
according to the eigenvalues of the monodromy, as well as between the orbits which differ 
by the discrete invariant (see Corollary |7.8| ). While all elliptic orbits with arbitrary values 
of the discrete invariant can be represented by Hill's operators with constant coefficients, 
just one of the hyperbolic and one of the Jordan block classes has such representatives, 
while others do not. 

Note that the bihamiltonian equations corresponding to elliptic orbits with different 
discrete invariants are almost the same: these are the CH equations with different co- 
efficients. It would be very interesting to see whether an analogous similarity holds for 
the hyperbolic and Jordan block orbits with different discrete invariants. 



6.6 Remark. When the symplectic leaves of the bracket {.,.}a are of codimension 
1, then the choice of a Casimir function is essentially unique for every A. (Any two 
Casimir functions for every fixed A are functionally dependent.) Therefore, the choice 
of the Poisson pair itself defines the bihamiltonian system (modulo the mentioned func- 
tional dependence of the initial Hamiltonian), provided that the symplectic leaves are 
hypersurfaces. 

This is indeed the case for the Virasoro coadjoint orbits discussed above, which 
have codimension 1 for a fixed cocentral value a. It turns out that a natural Casimir 
function hx{u{x){dxy,a) corresponding to the KdV Poisson pair on the dual space to 
the Virasoro algebra is the trace of the monodromy operator associated to Hill's operator 
—adl + u{x) — A^. It generates the first integrals of the KdV equation, see Remark |7.11 
in Appendix. Similarly, one can expand Casimir functions for the two other integrable 
cases, the CH and HS equations. 

Note that for orbits of higher codimension one can start with several Casimirs and 
consider several Lenard schemes to generate the sequences of Hamiltonians. 

6.7 Remark. A generic Virasoro coadjoint orbit Diff (S'^)/S'^ can be equipped with a 
complex structure and a two-parameter family of compatible (pseudo-) Kahler metrics 
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| Ki^ . 

This family of Kahler metrics has a simple origin: a generic Virasoro orbit has codi- 
mension 2, i.e., it is locally included in a two-parameter family of orbits, each equipped 
with its own symplectic structure compatible with the complex structure. Alternatively, 
one could consider a two-parameter family of symplectic structures on the same orbit, 
given by the Hamiltonian operators aSi^ + hdx- 

It turns out that the restriction of the two-parameter family of iJ^ ^ -metrics on vir* 
to a generic Virasoro orbit Diff (S"^) / 5*^ coincides with the family of Kahler metrics on the 
orbits. Proof is achieved by comparison with formula (7) in ||Kir2|| for those homogeneous 
metrics at one point of the orbit. 

This is yet another fact manifesting a special role of the -ff^^ -metrics in Virasoro 
geometry. 



7 Appendix: Classification of Virasoro orbits 

In this section we recall the classification result for Virasoro coadjoint orbits (see, e.g. 



Kir||, ||Seg|| or the book ||GK||). The dual spaces to the infinite-dimensional Lie algebras 



considered below are always understood as smooth duals, i.e. identified with appropriate 
spaces of smooth functions. 

A. Classification of quadratic differentials. We start with the non-extended group 
of difFeomorphisms of the circle. Let Diff(S'^) be the group of all orientation-preserving 
diffeomorphisms of 5*^ and let vect{S^) be its Lie algebra. 

Proposition 7.1 ||Kir|| The dual space vect{S^)* is naturally identified with the space of 
quadratic differentials {u{x){dx)'^} on the circle. The pairing is given by the formula: 

{u{x){dxY , v{x)dx)) = / u{x)v{x) dx 

for any vector field v{x)dx € vect{S^). The coadjoint action coincides with the action 
of a diffeomorphism on the quadratic differential: for a diffeomorphism (f G Diff (5) the 
action is 

Ad; : u {dxf ^ uiif) ■ ifl {dxf = u{ip) ■ {diff. 

Hence, for instance, if u{x) > for all x G S*^ then the square root ^yu{x){dx)^ 
transforms under a diffeomorphism as a differential 1-form. In particular, ^{u{x){dxY) = 
dx is a Casimir function (i.e., an invariant of the coadjoint action). One can 
see that there is only one Casimir function in this case, since the corresponding orbit 
has codimension 1 in the dual space vect{S^y. Indeed, a diffeomorphism action sends 
the quadratic differential u{x){dx)'^ to the constant quadratic differential C^dx)"^, where 
the constant C is the average value of the 1-form ^^u{x) dx on the circle: 

C = — / a/ u{x) dx. 
27r J51 
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On the other hand, if a differential u{x){dxy changes sign on the circle, then the 
integral j'' ^y\u{x)\ dx, evaluated between any two consecutive zeros a and b of the func- 
tion u{x), is invariant. In particular, the coadjoint orbit of such a differential u{x){dx)'^ 
has necessarily codimension higher than 1. 

7.2 Remark. In our study of the KdV equation we pick the "freezing" point in the dual 
space vect{S^)* to be C = 1/2. (Actually, we consider the dual space to the Virasoro 
algebra, but we choose the cocentral term equal to zero, so that the "freezing" point 
{uo{dxY, oq) = {{dxy/2, 0) belongs to the dual space to the Lie algebra of vector fields.) 
Other values of C give equivalent equations, different from the KdV by scaling only. 



B. Virasoro dual and Hill's operators. Let vir be the Virasoro algebra. We can think 
of its dual space as the space of pairs vir* = {{u{x){dxY,a)} consisting of a quadratic 
differential and a real number (cocentral term). It is more convenient, however, to regard 
such pairs as Hill's operators, i.e. differential operators —adl + u{x), as we will see below. 

Proposition 7.3 The Virasoro coadjoint group action is given by the formula 

Ad(^,b) : {u{dx)'^, a) i-^ {u{(p) ■ (pl{dx)'^ - aS{ip) {dxf, a) , (7.12) 

where 



2 



is the Schvuarzian derivative of(f. 



This group action on Hill's operators: 



Ad 



- adl + u{x) ^ -adl + u{^) ■ Lp\ - aS{Lp) 



has the following geometric interpretation (see, e.g., |peg|| , |[Kir|| , |Pvs|| ). Fix a = — 1 and 
consider Hill's operators of the form 9^ \u{x), x e . Let / and g be two independent 
solutions of the corresponding differential equation 

(a^ + M(x))y = o 



for an unknown function y. Although the equation has periodic coefficients, the solutions 
need not necessarily be periodic, but instead are defined over M. Consider the ratio 
r/ := //^ : M ^ MP\ 



Proposition 7.4 The potential u is (one half of ) the Schwarzian derivative of the ratio 
Tj: 

2 

Proof. Note that the Wronskian W{f,g) := fg^ — fxQ is constant, since it should 
satisfy the differential equation = 0. Here we normalize W by setting W = —1. This 
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additional condition allows one to find the potential u from the ratio rj. Indeed, first one 
reconstructs the solutions /, g from the ratio rj by differentiating: 



Therefore, 



Vx 



9 



fx9 - f9x -W 



9' 



1 



9' 



f = 9-V 



V 



/r]x ^"1]^ 

Given two solutions / and one immediately finds the corresponding differential equa- 
tion they satisfy by writing out the following 3 x 3-determinant: 



det 



y f 9 



Since / and g satisfy the equation y^rc + u ■ y = 0, one obtains from the determinant 
above that 

fx 9x 



u 



det 



The explicit formula for u expressed in terms of 7] turns out to be one half of the 
Schwarzian derivative of rj. □ 



Corollary 7.5 The Schwarzian derivative S{ri) is invariant with respect to a Mohius 
transformation rj [arj + b)/ [crj + d), where ad — be = 1. 

Proof. Indeed, for a given potential u the solutions / and g of the corresponding 
differential equation are not defined uniquely, but up to a transformation of the pair (/, g) 
by a matrix from 5*172 (M). Then the ratio rj changes by a Mobius transformation. □ 



Proposition 7.6 The Virasoro coadjoint action of a diffeomorphism (p on the potential 
u{x) gives rise to a diffeomorphism change of coordinate in the ratio rj: 

Proof. We look at the corresponding infinitesimal action on the solutions of the 
differential equation {d^ + u{x))y = 0. For (p{x) = x + ev{x) close to the identity, 
consider the infinitesimal Virasoro action of ip on the potential u{x): 

u ^ u + e ■ 5u, where 5u = 2uVx + u^v — -Vxxx, 

(cf. formula (|3.8|) for a = |). It is consistent with the following action on a solution y of 
the above differential equation: 

y v-^ y + f 5y, where 5y = -^yv^ + y^v. 
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The consistency means that {d^ + u + e ■ Su){y + e ■ 5y) = + (9(e^). 

Note that the action e ■ 5y = e ■ {—^yVx + yxv) is an infinitesimal version of the 
following action of the diffeomorphism (p{x) = x + ev{x) on y: 

ip: y t-^ y{ip){ip,^y^/'^. 

Thus solutions to Hill's equation transform as forms of degree —1/2. Therefore the ratio 
1] of two solutions transforms as a function under a diffeomorphism action. □ 

In short, to calculate the coadjoint action on the potential u one can first pass from 
this potential to the ratio of two solutions, then change the variable in the ratio, and 
finally take the Schwarzian derivative of the new ratio to reconstruct the new potential 

All of the above considerations of Hill's operators were local in x. To describe the 
Virasoro orbits, we now recall that u(x) is defined on a circle. 



Theorem 7.7 [^egj] , [[Kir|| The coadjoint Virasoro orbits (for a given cocentral term 
a ^ 0) are enumerated by the conjugacy classes in (S'L2(M) \ {id})/Z2, the universal 
covering of SL2{M.) without the identity and modulo the 'L2-action. 

Proof. For a periodic potential u the solutions of (9^ + u)y = are quasiperiodic. 
In other words, the boundary values of the fundamental set of solutions F := {f,g) on 
[0,27r] are related by a holonomy matrix M e SL2{R): F(27r) = F{0)M. (Similarly, 
one can consider the "projective solution" O := {ri,rix), which consists of the solution 
ratio r) and its derivative rj^ with the correponding holonomy Ai now in P S L2{M.) .) This 
holonomy matrix M changes to a conjugate matrix if xq = is replaced by any point 
xq G or if F is replaced by another system of solutions. 

Now regard the ratio r] = f/g as a map t] : [0,27r] — >■ MP^ describing a motion 
("rotation") along the circle MP^ ^ S^. One can see that the condition 77^. 7^ is 
equivalent to the condition W ^ on the Wronskian. Choosing the negative sign of 
the Wronskian, 1^ < 0, we can assume that the rotation always goes in the positive 
direction: rjx = —W/g"^ > 0. 

By a diffeomorphism change of the parameter x ^ <^{x), one can always turn the 
map rj : [0, 2-k] MP^ into a uniform rotation along MP^, while keeping the boundary 
values of r]{x) on the segment [0,27r] satisfying the holonomy relation B(27r) = O(0)A^. 
Furthermore, the number of rotations (the "winding number") for the map rj : [0, 27r] — > 
MP^ does not change under a reparametrization by a diffeomorphism Thus the orbits 
of the maps rj (or, equivalently, of the potentials {u{x)}) are described by the conjugacy 
classes of matrices in the universal covering of 5*172 (M). The choice in the sign of the 
Wronskian reflects the Z2-action on this universal covering. ^ 

Finally, note that the identity matrix in the universal covering 5*^2 (M) (or in its 
projectivization) cannot be obtained as a holonomy matrix for the maps rj : [0, 2n] — »• 
RP\ Indeed, any map r] starts at the identity and goes in the positive direction. Thus, 
no matter how slow the rotation, one always moves out from the identity. □ 
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Corollary 7.8 The Virasoro orbits in the hyperplane {—ad"^ + u{x) \ a = oq} C vir* 
with fixed a ^ are classified by the Jordan normal form of matrices in SL2{^) and 
a positive integer parameter, winding number. Matrices in the group SL2{M.) split into 



three classes, whose normal forms are 


the 


exponents of the following three classes: 





" /i " 
-/i 


, ii) 


' 



±1 " 



, and Hi) 


'00" 






in the corresponding Lie algebra 5/2(1^), see Fig. 4- The Virasoro orbit containing Hill's 
operator + u{x) has the codimension in the hyperplane {a = ao} C vir* that is 



equal to the codimension in 5*1/2 (M) of (the conjugacy class of) the holonomy matrix M 
corresponding to this Hill's operator. This codimension is 1 for the classes i) and ii), 
and it is 3 for Hi) independ of the integer parameter. 

7.9 Remark. Note that the class of real matrices whose exponents have the normal 
form i) decomposes into rotation matrices (with /i G iM.) and hyperbolic rotations (with 
fi E M.), the elliptic and hyperbolic cases, cf. Remark |6]^. In the paper we regard these 
cases as belonging to the same general class (and similarly we do not distinguish between 
the cases ±1 in type ii)), since we are interested in the algebraic (rather than geometric) 
question of constructing the corresponding integrable equations. 

One should also mention that the equality of codimension of the Virasoro coadjoint 
orbits in vir* and (the conjugacy class of) the corresponding holonomy matrices in in 
5*^2 (M) can be seen by checking the smooth dependence on a parameter in the above 
classification. (The versal deformations of the orbits can be given in terms of the Jordan- 
Arnold normal forms of the holonomy matrices depending on a parameter, cf. | |UK2| |.) 
Alternatively, one can find the dimension of the corresponding stabilizers, see ||Kir|| 



Regarded as homogeneous spaces, the orbits of type i) are often denoted by Diff (S*^) / S^, 
the notation Diff (S'^)/M^ stands for ii) (and sometimes for the case e iM in i)), and 
Diff(5'"'^)/S'L2(]R) corresponds to Hi). 




Figure 4: The points Mi, M2 and in shCR) (which is a local picture of SL2{M.)) correspond to the 
Virasoro orbits of the types i), ii) and Hi), respectively. 

7.10 Remark. For applications to bihamiltonian systems we would like to describe all 
points in vir* belonging to orbits of codimension at most 2. As we have shown above, 
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in the smooth dual there are no orbits of codimension or 1, as a is a Casimir function, 
and in each hyperplane the orbits are of codimension at least l.Q 

For ao 7^ the orbits are represented by the Hill's operators, whose holonomy matrices 
were classified above, while for = they are quadratic differentials. Our choices of 
representatives for the orbits of codimension 1 for a fixed qq (i-e. of total codimension 2 
in vir*) will be as follows. 

a) For a generic point representing the case i) above we take Hill's operator —d'^ + 
1/2 [uq = {dxY/2,aQ = 1). It corresponds to the differential equation {d^ — 1/2)?/ = 
and has the holonomy matrix diag (exp(7r-\/2), exp(— 7r-\/2)), the exponent of the type i). 
This freezing point corresponds to the CH equation. 

b) The matrix of type ii) can be encountered in a generic 1-parameter family of 
matrices in 3/2 (K)- Its exponent, a Jordan 2 x 2-block with the eigenvalue 1, can be 
represented as the holonomy matrix by the Hill's operator —d"^ [uq = 0,ao = 1). The 
latter point in vir* corresponds (after an appropriate reduction) to the HS equation. 

c) The hyperplane ao = in the Virasoro dual is the dual space vect* to the non- 
extended Lie algebra of vector fields. The orbits of codimension 1 in the space vect* are 
represented, e.g., by the quadratic differential ^{dxY (i.e., by the point (ao = 0,mo = 
{dxY/2)) as we discussed in Section A. Freezing the Poisson structure at the latter point 
leads to the KdV equation. 

The above three cases are described in Corollary 



7.11 Remark. Recall that the holonomy matrix M of Hill's operator d"^ + u{x) 
changes to a conjugate one under the Virasoro action. This implies that h{dl + u{x)) = 
log(trace M) is a Casimir function on vir*. One can use it to generate the KdV hierarchy 
via the Lenard scheme described in Theorem |6.3| . 

Recall that for the KdV equation the freezing point for the constant Poisson structure 
{., .}o is (ao = 0,Mo = {dxY/2). Therefore, the Casimir function of the bracket {., .}a : = 
{■) - jip — .}o has the form hx{dl+u{x)) := log(traceMA), where Mx is the holonomy 
of the operator d"^ + u{x) — X^. The expansion of the function hx in A produces the first 
integrals of the KdV equation: 



hx ^ 2nX - ^ c„/i2„-iA^ 



2n 
n=l 



where 



hi = I u{x)dx, hs = I u^{x)dx, h^ = I [ u^{x) ~ -{ux{x))^ ] dx, 
'51 Js^ Js^ 



2 



and Ci = 1/2, c„ = {2n — 3)!!/(2"r;,!) for n > 1. One can see that the Hamiltonian 
/13 is quadratic in u and coincides with the "energy" Hamiltonian of the KdV equation, 
regarded as an Euler equation. (Note that the KdV Hamiltonians hj are differential poly- 
nomials whose degree increases with j. The latter follows from the recurrence relation 

^There exist Virasoro orbits of codimension 1 if in the dual space vir* besides smooth elements we 



also admit singular ones, cf. |Wit]. In this paper we consider the classification of the smooth dual 
elements only. 



24 



{^2j+i5 /}o + {^2j-i) /}lp = for Hamiltonians hj (cf. Theorem |6.3|) for the constant 



and hnear Poisson brackets on vir* .) Much more details on the KdV structures can be 
found in |PZ| . 

Similar computations can be done for CH and HS, the other two equations considered 
in this paper, cf. ||BSS|| . 



Acknowledgements 

We are grateful for the hospitality to the Isaac Newton Institute (Cambridge, UK) and 
the Fields Institute (Toronto, Canada), where most of this work was done. We are 
indebted to P. Pushkar and D. Novikov for drawing the figures. Special thanks go 
to I. Zakharevich for many useful remarks, as well as for illuminating discussions on 
Theorem ff?7. 

The work of B.K. was partially supported by an Alfred P. Sloan Research Fellowship 
and by an NSERC research grant. The work of CM. was supported in part by NSF 
Grant DMS-9970857. 



References 

[Arn] V. Arnold, Sur la geometrie differentielle des groupes de Lie de dimension infinie et 
ses applications a I'hydrodynamique des fluides parfaits, Ann. Inst. Fourier (Greno- 
ble) 16 (1966), 319-361; Mathematical Methods of Classical Mechanics, Springer- 
Verlag, New York 1989, xvi+508 pp. 

[AK] V. Arnold and B. Khesin, Topological Methods in Hydrodynamics, Springer- Verlag, 
New York 1998, xv+376 pp. 

[BSS] R. Beals, D.H. Sattinger and J. Szmigielski, Acoustic scattering and the extended 
Korteweg-de Vries hierarchy. Adv. in Math., 140 (1998), 190-206; Inverse scatter- 
ing solutions of the Hunter-Saxton equations, preprint (2001). 

[Bre] Y. Brenier, Minimal geodesies on groups of volume-preserving maps and generalized 
solutions of the Euler equations. Comm. Pure Appl. Math. 52 (1999), no. 4, 411- 
452. 

[CH] R. Camassa and D. Holm, An integrable shallow water equation with peaked solu- 
tions, Phys. Lett. Rev. 71 (1993), 1661-1664. 

[FF] A. Fokas and B. Fuchssteiner, Symplectic structures, their Bdklund transformations 
and hereditary symmetries, Phys. D 4 (1981/82), 47-66. 

[GZ] I.M. Gelfand and I. Zakharevich, Spectral theory of a pencil of third-order skew- 
symmetric differential operators on S^, Funct. Anal. Appl. 23 (1989), 85-93. 

[GR] L. Guieu and C. Roger, The Virasoro algebra and Virasoro group, monograph in 
preparation. 



25 



[GS] V. Guillemin and S. Sternberg, On collective complete integrability according to the 
method of Thimm, Ergodic Theory Dynam. Systems 3 (1983), no. 2, 219-230. 

[HS] J.K. Hunter and R. Saxton, Dynamics of director fields, SIAM J. Appl. Math. 51 
(1991), no. 6, 1498-1521. 

[HZh] J. Hunter and Y. Zheng, On a completely integrable nonlinear hyperbolic varia- 
tional equation, Phys. D 79 (1994), 361-386. 

[Kir] A. Kirillov, The orbits of the group of diffeomorphisms of the circle, and local 
Lie superalgebras, Funct. AnaL AppL 15 (1981), 135-136; The orbit method, II: 
Infinite- dimensional Lie groups and Lie algebras, Contemp. Math, 145 (1993), 
33-63. 

[Kir2] A. Kirillov, Kahler structure on the K -orbits of a group of diffeomorphisms of the 
circle., Funct. Anal. Appl. 21:2 (1987), 42-45. 

[Kru] M. Kruskal, Nonlinear wave equations, in "Dynamical Systems, Theory and Ap- 
plications", Springer- Verlag, Heidelberg 1975. 

[Ma] S. Manakov, A note on the integration of Euler's equations of the dynamics of an 
n-dimensional rigid body, Funct. Anal. Appl. 10 (1976), 93-94. 

[Mi] G. Misiolek, A shallow water equation as a geodesic flow on the Bott-Virasoro 
group, J. Geom. Phys. 24 (1998), 203-208; Classical solutions of the periodic 
Camassa-Holm equation, to appear in Geom. and Funct. Anal. 

[Ovs] V. Ovsienko, On the Denogardus great number and Hooke's law, Kvant (1989), 
8, 8-16; Kvant selecta: Algebra and Analysis, H, Math. World, 15, Amer. Math. 
Soc, Providence, RI, (1999) 153-159. 

[OK] V. Ovsienko and B. Khesin, The (super) KdV equation as an Euler equation, Funct. 
Anal. Appl. 21:4 (1987), 81-82. 

[0K2] V. Ovsienko and B. Khesin, Symplectic leaves of the Gelfand-Dikii brackets and 
homotopy classes of non- degenerate curves, Funct. Anal. Appl., 24:1 (1990), 33-40. 

[Seg] G. Segal, Unitary representations of some infinite- dimensional groups. Comm. 
Math. Phys. 80 (1981), 301-342; The geometry of the KdV equation, in "Topologi- 
cal Methods in Quantum Field Theory" (Trieste, 1990), Internat. J. Modern Phys. 
A 6 (1991), 16, 2859-2869. 

[Th] A. Thimm, Integrable geodesic flows on homogeneous spaces, Ergodic Theory Dy- 
nam. Systems 1 (1981), 491-517. 

[Wit] E. Witten, Coadjoint orbits of the Virasoro group. Comm. Math. Phys. 114 (1988), 
1-53. 

[Za] I. Zakharevich, manuscript in preparation. 



26 



